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Abstract 



We review recent results concerning finite size corrections to the Ising 
model free energy on lattices with non-trivial topology and curvature. From 
conformal field theory considerations two distinct universal terms are ex- 
pected, a logarithmic term determined by the system curvature and a scale 
invariant term determined by the system shape and topology. Both terms 
have been observed numerically, using the Kasteleyn Pfaffian method, for lat- 
tices with topologies ranging from the sphere to that of a genus two surface. 
The constant term is shown to be expressible in terms of Riemann theta func- 
tions while the logarithmic correction reproduces the theoretical prediction by 
Cardy and Peschel for singular metrics. 
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I. INTRODUCTION 



Finite size corrections provide the shortest bridge between the critical behavior of a 
lattice model and the conformal field theory description of its universality class. The simplest 
example of such a connection being the dependence of finite size corrections on the model 
central charge for a strip geometry^. In more general geometries, the universal part of the 
free energy will depend sensitively on the lattice shape, topology and curvature. For a two 
dimensional system, characterized by a length scale L, the free energy at criticality and for 
fixed shape has a large L expansion^ of the form 

F = /oL2 + /;,L + ClnL + L' + o(l). (1.1) 

where C and D are universal, the logarithmic term coefficient C being determined by the 
system curvature and the constant term D by its topology and shape. 

In a series of recent papers^^^ the author and collaborator have studied these terms for 
the Ising model on lattices with topologies ranging from the sphere to that of a genus two 
surface. In this talk we will review these results using as examples the square lattice with 
spherical topology. Fig. 1, and the square lattice embedded on a genus two surface. Fig. 2. 

While it would be desirable to explicitly evaluate C and D in the lattices scahng hmit, 
progress along this direction has been hindered by difficulties in evaluating the free energy. 
Closed form expressions for the Ising model free energy are known only for a small number 
of finite size lattices: the torus^, the cylinder^ and some non-orientable^ geometries, where 
translational invariance permits the use of Fourier transforms. Progress can be made however 
by using the Kasteleyn formalism to express the partition function as a linear combination of 
determinants of adjacency matrices. These determinants can then be evaluated numerically 
for large lattice sizes and arbitrary lattice geometry. 

The Ising model on lattices with the topology of the sphere has been studied 
previously^"^^ and the coefficient C was obtained for a restricted class of hexagonal lattices^^. 
Topology dependent terms are also observed, for spherical topology, in the context of 
percolation^^. 

The paper is organized in the following way: 

In Section II we review the Kasteleyn method for the evaluation of the Ising model 
partition function on general two dimensional lattices. 

In Section III we discuss the constant term D for lattices embedded on a genus g sur- 
face. Careful numerical studies showed that D can be expressed in terms of Riemann theta 
functions, defined below in the text, with half-integer characteristic k 

Dc^b Y.\e[k]m)\- (1-2) 
fc=i 

This result is reminiscent of multi-loop calculations in string theory^^'^^ with D reproducing 
the modular invariant partition function of the c = | conformal field theory on a genus 
g surface. The lattice shape and couphng constants determine the g x g matrix Q. This 
matrix is the period matrix of a Riemann surface related with the lattice continuum limit. 
In Section IV we study the dependence of Q on the lattice shape using a discrete formulation 
of Riemann surface theory. 
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In Section V we discuss the logarithmic term C for both spherical and higher genus 
lattices. From a conformal field theory point of view it is known^ that this term is sensitive 
to the smoothness of the system. For a smooth system it is related with the total curvature 

by 

C = -^cx (1.3) 

with X = '2—'2g being the system Euler characteristic and c the central charge. By opposition, 
the contribution from a conical singularity with deficit angle e = 27r — ^ is 

Ce = ^{l-{2n/9r) (1-4) 

and not a term proportional to e as one would expect from a delta function in the curvature. 
For the scaling limit of a specific lattice the distinction between the two cases depends on the 
order in which the limit of singular curvature and the thermodynamic limit are taken. We 
will show that, for the lattices of Figs. 1 and Fig. 2, C is given by the sum over all conical 
singularities in the lattice of the contribution (1.4). Similar results have been reported by 
Gonzalez^^ . 

Finally in Section VI we present our conclusions and discuss open problems and avenues 
for future research. 



II. THE KASTELEYN FORMALISM 

Consider a lattice, or more properly a graph, G drawn without superposition of edges on 
a surface of genus g. Assigning a coupling constant to each edge and placing an Ising spin 
at each vertex we can defined the Ising model on G. Kasteleyn^^'^^ showed that its partition 
function can be expressed as a linear combination of the Pfaffians of 4^ antisymmetric 
matrices. This result provides the starting point for the study of the Ising model on non- 
trivial geometries. In this section we provide a brief description of the Kasteleyn method, 
for further details and proofs we refer the interested reader to references^'^^'^^. 

Given a graph embedded on a genus g surface, the Kasteleyn method to obtain the Ising 
model partition function consists of essentially three steps: 

First step: replace each vertex of the original lattice by a cluster of vertices, called the 
decoration graph, that depends on the vertex coordination number. A possible choice of 
the decoration graph for vertices with three, four and five nearest neighbors is shown in 
Fig. 3. The order by which the exterior edges connect to the decoration is not important, 
meaning that the decoration graphs can be rotated. The resulting decorated lattice will be 
denoted by Gd- Its edges are assigned weights: 1 to the internal edges of the decorations 
and a weight w — tanhX to the edges inherited from the lattice G, with K being the Ising 
coupling constant on that edge in units of KbT. 

Second step: orientate the edges of Gd by assigning to each edge a direction, represented 
graphically by an arrow. The internal edges of the decorations have already an orientation 
as given in Fig. 3. The remaining edges, inherited from the original lattice G, are oriented 
according to the Kasteleyn rule: in such a way that all lattice faces have an odd number of 
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clockwise oriented edges. An example of an orientation satisfying the Kasteleyn condition 
is shown in Fig. 4 for a genus two lattice. 

For any lattice there are many different edge orientations that satisfying the Kasteleyn 
condition. Two lattice edge orientations are said to be equivalent if they can be related one 
to the other by a series of arrow reversals, in which the orientations of all the edges meeting 
at a given vertex are reversed. It can be shown^'-'^^'-'^^ that for a genus g lattice there are 
precisely 4^ un-equivalent Kasteleyn orientations, matching the number of spin structures 
on the corresponding continuum free fermion description. Let ai,bi with i = l,...,g be 
closed paths on the dual lattice forming a canonical basis of the embedding surface first 
homology group; then the 4^ un-equivalent Kasteleyn orientations can be generated from an 
initial Kasteleyn orientation by reversing the orientation of the edges crossed by a choice of 
such loops. 

Third step: label the vertices of the decorated graph Gd with an integer from 1 to A^. 
To each Kasteleyn edge orientation associate a M. x A4 adjacency matrix A with entries 
Aij that vanish if vertices i and j are not connected by an edge and take a value ±z if 
vertices i and j are connected by an edge of weight z, the sign being determined by the edge 
orientation. Schematically 



Aij — < 



z Hi — ^ j 

-z iij^i . (2.1) 
otherwise 

The Ising model partition function on the graph G is then given by 

Z{K) = 2^- (cosh X)^- - ^a,Pf ^,(X) (2.2) 

^ i=i 

were Ny and Ne are the number of vertices and edges in G and we assume that all edges have 
the same coupling constant K. The sum runs over representatives of the 4^ un-equivalent 
Kasteleyn orientations. The take values ±1 and are completely determined by the arrow 
parity of non-trivial topology loops a^, bi along lattice edges, see'^ for details. 

This completes our description of the Kasteleyn method. We have reduced the evaluation 
of the Ising model partition function to the calculation of determinants of adjacency matrices, 
using that Pf A^ — ^/detAi for an antisymmetric matrix. 



III. THE CONSTANT TERM D ON HIGHER GENUS LATTICES 

Consider the genus two lattice shown in Fig. 2, it can be seen as a torus with an additional 
handle on the bulk and it is characterized by five integer sizes Mj. Locally it is equivalent to 
the flat square lattice except around the two octagonal faces. These faces in the continuum 
limit will correspond to conical singularities on a otherwise flat system. To study the flnite 
size corrections we evaluate the free energy on a sequence of lattices with flxed shape and 
increasing size. This is done by taking lattice dimensions of the form Mj = L with fixed 
rrii and increasing L. The coupling constants are fixed for all edges to the square lattice 
isotropic critical value 
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smh2Kc = 1. 



(3.1) 



Following the discussion of the previous section, the partition function on a genus two 
lattice can be evaluated in terms of sixteen determinants of adjacency matrices. These 
matrices and the associated Kasteleyn edge orientations can be labeled as A{na^ , n^^ . n^,^ , J 
with Tlx = 0, 1. The starting edge orientation A(OOOO) is shown in Fig. 4 and an orientation 



with n^: = 1 is obtained from the corresponding orientation with 



by reversing 



the orientation of all the edges crossed by the cycle x. To make connection with the theta 
function characteristics and allow for more compact equations we will also use the alternative 
notations 



A, 



(3.2) 



with the integer label given by i = 16 — Sn^^ — 4n^^ — 2na2 

In terms of these orientations, the Ising model partition function on the genus two lattice 
is given by 



T' 



Pi -P2-P3-PA- 

+Pio + Pn + Pi2 



P5 + Pe + P7 + P8- P9 

- Pi3 + Pu + Pi5 + Pie) 



(3.3) 



where Pi = Pf Ai and ao is the size dependent pre-factor in (2.2). Closed form expressions 
for these Pfaffians are difficult to obtain, since the lattice is not translationally invariant, 
but they can be evaluated numerically for relatively large lattice sizes. 

To factor out the bulk term /o and the logarithmic correction C we can, as in the 
toroidal case^'^, consider ratios of determinants of adjacency matrices to the largest among 
the sixteen. These ratios are found to converge to well defined values for large lattice size 
Ny, see Fig. 5 for an example. More precisely, we find that in the L — > 00 hmit the ratios 
of determinants satisfy at criticality 
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(3.4) 



for the 16 combinations of c^; = 0, 1 and with the label i defined below equation (3.2). The 
genus two Riemann theta functions are defined by the quickly converging series 



e 



a 
/3 



(z,fi) 



^ exp 

neZ2 



m{n + OLfQ.{n + a) + 27ri(n + a)^(z + (3) 



(3.5) 



where a, /3, z and n are 2-vectors half-integers, complex numbers and integers respectively. 
Six of the sixteen theta functions in (3.4) are odd functions of z and therefore vanish at the 
origin. The corresponding ratios of determinants are found to vanish in the large Ny limit. 

The 2x2 period matrix f2 is a symmetric, complex valued matrix with a positive definite 
imaginary part. It is determined by the lattice shape and coupling constants, its value in 
(3.4) can be obtained directly by numerically fitting the ratios of determinants or from first 



5 



principles, as we will see in the next section, by using discrete holomorphy methods. For 
locally square lattices, fl is found to be purely imaginary 



fill Qi2 
^12 ^22 



(3.6) 



this is not the case for triangular lattices where the period matrix is in general complex^. The 
same property is observed on toroidal lattices where the modular parameter is imaginary 
for squared lattices^'^ and complex for triangular lattices^^ . 

In Table I these results are illustrated for the lattice with shape (rrii) = (1, 1, 1, 1, 1). The 
fifteen determinant ratios are shown, both the largest lattice size evaluation and the L ^ oo 
extrapolation obtained using a fit with a quadratic polynomial in l/Ny. These ratios are 
compared with ratios of theta functions 



©(16-8di-4ci-2d2-C2)(^) — ® 



ci/2 C2/2 
di/2 d2/2 



(o,n)/e 1 1 (o,n) 



(3.7) 



with a period matrix obtained by a suitable numerical fitting procedure. The precision to 
which the two sets of numbers agree is remarkable, typically a precision from 10~^ to 10~^. 



IV. DISCRETE HOLOMORPHY AND THE PERIOD MATRIX 

The matrix Q in Eq. (3.4) can be seen as the period matrix of a Ricmann surface related 
with the lattice continuum limit. In this section we study the dependence of this matrix on 
the lattice shape using a discrete formulation of holomorphy. The basic idea is to formulate 
Riemann surface theory on a discrete setting by using finite difference operators^'^'^^'^^. 
These operators act on quantities defined on the lattice p-elements: vertices, oriented edges 
and faces that we will call respectively the lattice functions, differentials and volume forms. 

A lattice function / is determined by its value on the lattice vertices and can be rep- 
resented by a A''y- vector f[n] : n = 1, . . . ,Nv after an integer labeling of the lattice Ny 
vertices is chosen. 

A lattice differential w is determined by its value on the lattice oriented edges. Referring 

to a fixed drawing of the lattice, wc define horizontal edges to be oriented from left to 
right and vertical edges from bottom to top. A lattice differential is then represented by a 
Ny X 2-matrix w[n|p] : n = 1, . . . , Ny] p = 1,2 where [n|l] stands for the horizontal edge 
immediately right of vertex n and [n|2] for the vertical edge immediately above n, see Fig. 
7. The integral of a lattice differential w along a path of lattice edges C is the sum of the 
values that w takes on the edges included in C 

w = ^ ±w[n\p] (4.1) 

^ [n|p]GC 

with a minus sign for edges with opposite orientation to that of the path. 

A lattice volume form rj takes values on the lattice faces and is represented by a {Ny — 2)- 
vector ri[q] : q — 1, . . . , Ny — 2 after a labeling of the Ny — 2 lattice faces is chosen. The 
integral of a lattice volume form rj over a given lattice area A is the sum of the values that 
Tj takes on all faces included in A. 
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It is important to relate the labeling of vertices, edges and faces. Vertices and edges are 
already related by the labeling introduced above, to relate vertices with faces we introduce 
some additional notation: if g is a lattice squared face then qi is its lower left vertex, if q is 
an octagonal face then qi and q2 correspond to the two vertices that can be seen as its lower 
left vertices, see Fig. 2. Conversely for a given vertex n we denote by h the lattice face of 
which n is a lower left vertex. 

We now define finite difference operators acting on these quantities, which are the lattice 
versions of the exterior derivative, the co-derivative and the Hodge star operator, see Fig. 
6. 

The lattice exterior derivative d is a linear operator defined by 

(c^/)[n|l] = /[right(n)]-f[n] (4.2) 
{df) H2]=/[up(n)]-f[n] 
{dw) [q] = X: {wim + «;[right(qi)|2] - w[up(qi)|l] - w[qi|2]) 

i 

with i=l,2 for octagonal faces and i=l for squared faces. The functions right(n), left(n), up(n) 
and down(n) give the label of the vertex immediately right, left, above and below of the vertex 
n. 

The lattice co-derivative S is the operator defined as 

(Sw) [n] = w[left(n)|l] - w[n|l] + w[down(n)|2] - w[n|2] (4.3) 
{S rj ) [n|l] = ri[h] — ?7[down(n)] 
{St]) [n\2] = ?7[ieft(n)] - 7][n] 

(4.4) 

where left(n) stands for the face given by the tilde of the vertex left(n) and similarly for 
down(n). 

The discrete Hodge star acting on lattice differentials is defined as 

(^)[n|l] = -w[down(n)|2] (4.5) 
(^)[n|2]= 'y;[left(n)|l]. 

These discrete operators are defined in such a way that they satisfy most of the usual 
properties of their continuum counterparts: the lattice exterior derivative satisfies a discrete 
version of Stokes theorem; we have that dd — SS — and the two operators are adjoint 
of each other under the trivial inner product of lattice functions, differentials and volume 
forms. These properties are exact for all lattices at finite size. Such is not the case for some 
important properties of the discrete Hodge star operator that, as we will see below, are only 
satisfied in the limit of large lattice size. 

In exact analogy with the continuum definitions, a lattice differential is said to be har- 
monic if it is both closed and co-closed, 

w is harmonic = i idw)[q] = 0, q = 1, . . . , Np 

w is narmomc - j ^ q, n = 1, . . . , ' ^^"^^ 
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where the conditions must be met for all Np lattice faces and all Ny lattice vertices. In 
general, on a genus g lattice, the system of equations (4.6) has 2g linearly independent 
solutions. Notice that the number of of unknowns w[n|p] minus the number of equations is 
already 

NE-Nv-NF^-X^2g-2 (4.7) 

where Np is the number of lattice faces and x is the Euler characteristic of the lattice. 
Two additional solutions are provided by constant differentials with independent vertical 
and horizontal components. That these are the only solutions follows from the fact that a 
harmonic differential, discrete or continuum, is completely determined by its periods along 
the lattice, or surface, non-trivial loops and there are only 2g such loops in a genus g lattice. 

Finally a lattice differential is said to be holomorphic if it is harmonic and it satisfies 
icw — —iw. We then need to require the lattice Hodge star to be an endomorphism on the 
space of harmonic differentials and to satisfy ick = —1. Unfortunately these two properties 
are only approximately verified on a finite size lattice, see^ for details. Except for some 
specially symmetric lattices, one expects these properties to be exact only in the large Ny 
limit. 

As in continuum Riemann Surface theory, the lattice period matrix can be defined it 
terms of the lattice holomorphic differentials. We start by solving the linear system of 
equations (4.6) for the lattice harmonic differentials. There will be 2g independent solutions 
spanning a 2g linear space that can be decomposed, almost precisely, in two sub-spaces of 
holomorphic and anti-holomorphic lattice differentials by means of the projection operators 

P = (1 + i^) /2 and P = (1 - i^) /2. (4.8) 

We then look for a basis of the holomorphic sub-space {Ffe : k — 1 . . . ,g} with normalization 

/ ri^5ki (4.9) 

Jak 

[ ri = Qki (4.10) 

where the integrals are along a choice aj,bj of closed paths of lattice edges representing a 
basis of the embedding surface first homology group. 

Equation (4.10) gives our evaluation of the finite size period matrix. Detailed numerical 
studies^''^ showed the resulting matrix to be a complex, symmetric matrix with a positive 
definite imaginary part. For any lattice of the form shown in Fig. 2, the fixed shape Ny — > oo 
limit of this matrix reproduces with excellent precision the period matrix Q in equation (3.4). 
The lattice period matrix depends however on the choice of the first homology group basis 
but different period matrices for the same lattice will be related by a modular transformation 
on the limit of large lattice size. 

In Fig. 5 four ratios of determinants of adjacency matrices are plotted, in function of the 
lattice size A^^, together with theta function ratios (3.7) for period matrices evaluated with 
two different basis of the first homology group. The basis A is similar to the one shown in 
Fig. 4 while the basis B has the a.j and 6j loops interchanged from Fig. 4. We see that, for 
each ratio, the three values converge to a common value on the large Ny limit. In Table I 
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numerical values are given for the determinant and theta function ratios. The theta function 
ratios are shown both for a period matrix fitted numerically to the ratios of determinants 
and the L —>■ oo extrapolated lattice period matrix evaluated using the procedure described 
in this section. Numerical values for the entries of the two period matrices agree within 1%, 
the difference being due to the difficulty in the obtaining the L — > oo extrapolations. 

The method introduced in this section to evaluate lattice period matrices can be gener- 
alized to the anisotropic squared lattice^ and the general triangular lattice^. 



V. THE LOGARITHMIC TERM C 

We now consider the logarithmic term for both the spherical lattices of Fig. 1 and the 
genus two lattice of Fig. 2. The toroidal lattice will also be considered for illustration 
purposes. The Ising model free energy can be expressed in terms of the determinants of 
one, four and sixteen adjacency matrices for the spherical, toroidal and genus two lattice 
respectively. Again we will consider sequences of lattices with fixed shape rrii and increasing 
size Mi = rriiL at the isotropic critical coupling (3.1). For the higher genus lattices we 
learned from (3.4) that the different adjacency matrix determinants converge rapidly to a 
common bulk term times a topology and shape determined factor, 

Pf A(iVv) = e, Z{Nv) (5.1) 

where the ©j are constants for large enough lattice size Ny- It follows that since the 
logarithmic correction is entirely due to the Z term, 

-F = ln(^^£«,e,j+lnZ, (5.2) 

it can be equivalently evaluated from any non-vanishing Pfaffian. For this purpose it is 
convenient to introduce the auxiliary quantities 

HNv) - - In (^a.Pf A,) = F{Ny) + In ^) (5-3) 

where i labels one of the 4^ un-equivalent Kasteleyn orientations. In the thermodynamic 
limit Fi differs from F only by a constant. 

The logarithmic correction C can be found by subtracting to the lattice free energy the 
leading volume term /q, that is known^'' to be 

/o = 2G/7r+^ln2 (5.4) 

for the squared lattice, with G being the Catalan constant. In Fig. 8 the residual free 
energy F — f^Ny is plotted as function of the logarithm of the lattice size, (logA^y)/2, for 
lattices with the various geometries. A clear linear behavior is observed and the difference 
between positive, zero and negative curvature is patent for the spherical, toroidal and genus 
two lattice. 
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For each lattice we evaluate a series of values converging to C by doing linear regressions 
on sets of four lattices with consecutive increasing sizes. Examples of values obtained are 
shown in Table 11. We find convergence to a well defined value of C, the rate of convergence 
depending on the shape of the lattice, being fastest for the more symmetric lattices. The fact 
that the two spherical lattices, the cube and the L-shaped lattice, have different logarithmic 
corrections C already points out that this correction cannot be explained by the smooth 
metric contribution equation (1.3) but it is given by a sum over the conical singularities 
contributions of equation (1.4). The Cq corrections for vertices where three, five and eight 
squared faces meet are hsted on Table III. While the singularities of the genus two lattice 
are located at the octagonal faces we can by duality see them as vertices where eight faces 
meet. The total logarithmic correction C is obtained by summing the Cg contributions for 
all singularities occurring in a lattice. As shown in Table II, this reproduces the numerical 
values of C with high accuracy. The validity of equation (1.4) for all the lattices considered 
is clearly estabhshed from the numerical results. 

It should be remarked that the conical singularity behavior is related with the large scale 
structure of the lattice and not the small scale structure. In^ it was shown that a vertex 
where more, or less, than four faces meet does not necessarily constitute a conical singularity. 
Conical singularity behavior on the finite size corrections is observed when a regular lattice 
is folded at large scale into a conical shape spanning an angle 9. 

VI. CONCLUSIONS 

In this talk we reviewed recent research on finite size corrections, the universal terms, of 
the Ising model free energy in lattices with curvature. While these terms are well understood 
from a conformal field theory point of view they are not usually studied in the scaling limit 
of specific lattice geometries. Using the Kasteleyn method and numerical evaluation of the 
adjacency matrix determinants it was possible to study a number of different geometries 
consisting of regular lattices folded around conical singularities. 

The scale invariant term D was shown to be a modular invariant quantity expressed in 
term of Riemann theta functions, the dependence on the lattice shape and coupling constants 
being realized through a period matrix Vt. This result provides a lattice regularized picture 
of the classical winding part of the partition function in multi-loop calculations in string 
theory^^'^^. The study of the exact dependence of the period matrix on the lattice shape 
provides also a testing ground for ideas of discrete holomorphy. 

On the other hand the study of the logarithmic term C gave a first answer to the long 
standing question, posed by Cardy and Peschel in^, to which of the two forms (1.3) or (1.4) 
for the logarithmic correction is observed in specific lattice models. We found that for the 
regular lattices with conical singularities considered in this paper, the logarithmic correction 
is given by a sum over conical singularities contribution of the form (1.4). 

It is natural to ask what kind of lattice would have a logarithmic term given by equation 
(1.3), the smooth metric contribution. Such a lattice should have its curvature spread 
over the all lattice with the local curvature being non-zero everywhere but vanishing in 
the continuum limit. It is not clear how this can be accomplished or what the criticality 
condition should be for such a lattice. In this paper we have discussed a number of non- 
usual lattices boundary conditions but it seems that to obtain (1.3) one must also consider 
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lattices with non-usual structure or a non-usual limiting procedure. One thing is clear, since 
a continuum result must have a well defined regularization procedure, such lattices must 
exist. 

Acknowledgments 

This work was partially supported by the EU grant HPRN-CT-1999-00161. Most of the 
work described in this talk was done in collaboration with Prof. Barry McCoy. 



11 



REFERENCES 



1 H.W.J. B15te, J.L. Cardy and M.P. Nightingale, Phys. Rev. Lett. 56 (1986) 742. 

2 J.L. Cardy and L Peschel, NucL Phys. B300 (1988) 377. 

3R. Costa-Santos and B.M. McCoy, NucL Phys. B623 [FS] (2002) 439, hep-th/0109167 
^R. Costa-Santos and B.M. McCoy, J. Stat. Phys. 112 (2003) 889, cond-mat/02 10059. 
^ R. Costa-Santos, cond-mat/0306396. 
^B. Kaufman, Phys. Rev. 76, 1232 (1949). 

B. M. McCoy and T.T. Wu, The two-dimensional Ising model (Harvard Univ. Press, Cam- 
bridge, Massachusetts, 1973). 

«W.T. Lu and F.Y. Wu, Phys. Rev. E 63 (2001) 26107. 

9 J. Gonzalez, M.A. Martin-Delgado, preprint PUPT-1367, hep-th/9301057. 

°0. Diego, J. Gonzalez and J. Salas, J. Phys. A 27 (1994) 2965. 

^ Ch. Hoelbling and C.B. Lang, Phys. Rev. B 54 (1996) 3434. 

2 J. Gonzalez, Phys. Rev. E61 (2000) 3384. 

3 CD. Lorcnz and R.M. Ziff, Physica A 296 (2001) 1. 

^K. Kancda and Y. Okabe, Phys. Rev. Lett. 86 (2001) 2134. 
^W. Janke and R. Kenna, Phys Rev. B65 (2002) 064110. 

^L. Alvarez-Gaume, G. Moore and C. Vafa, Commun. Math. Phys. 106 (1986) 1. 
^E. Verlinde, H. Verlinde, NucL Phys. B288 (1987) 357. 

^ P.W. Kastclcyn in Graph Theory and Theoretical Physics, ed. F. Harary, 43 (Academic 

Press, New York, 1967). 
^N.P. Dolbihn, A.S. Mishchenko, M.A. Shtan'ko, M.I. Shtogrin and Y.M. Zinoviev, Func. 

Anal, and Applic. 30 (1996) 163. 
'°A.E. Ferdinand and M.E. Fisher, Phys. Rev. 185 (1969) 832. 

C. Nash and D. O'Connor, Phys. Rev. Lett. 76 (1996) 1196. 
'2 C. Mercat, Commun. Math. Phys. 218 (2001) 177. 

'^C. Mercat, math-ph/0111043. 



12 




FIG. 1. Two lattices with spherical topology: the cube and the L-shaped lattice, characterized 
by integer sizes Mj. Conical singularities are related with the corners where three or five faces 
meet. 
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FIG. 2. The genus two lattice can be seen has a toroidal lattice with an additional handle in 
the bulk. The boundary edge identifications are given by the letters. The two octagonal faces are 
marked in dashed line. 
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FIG. 4. The ^(0000) Kasteleyn edge orientation of the decorated genus two lattice with integer 
sizes (Mi) = (2, 4, 3, 3, 2). The remaining fifteen Kasteleyn orientations are obtained from this one 
by reversing the orientations of all edges crossed by a selection of the ai,bi loops. 
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FIG. 5. Ratios of determinants for the genus two lattice with shape (mj) = (11111) in 
function of the number of lattice points N and the theta function ratios 9f for period matrices 
evaluated using two different first homology group basis, A and B. The ratio -R13 corresponds to 
an odd characteristic theta function and vanishes in the large Ny limit. 
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FIG. 6. The lattice functions /, differentials w and volume forms 7] and the finite difference 
operators acting on them. 
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FIG. 7. The labeling of vertices, edges and faces for squared and octagonal faces. 
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FIG. 8. The residual free energy in function of the logarithm of the lattice size Ny for var- 
ious lattices. For the torus and genus two lattice only one Pfaffian term Fi is considered. The 
uncertainties on the values are smaller than the symbols used. 
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TABLES 



i 


it^(4A'y = 5184) 


Til /'AT \ 

Ri\Ny oo) 




^»(eval) 


1 


0.066121 


0.066386 


f\ francs 

0.066378 


0.066794 


2 


0.000012 


1.4 X 10 ° 


0.000000 


0.000000 


3 


0.000009 


1.3 X 10~^ 


0.000000 


0.000000 


4 


0.000009 


1.3 X 10^^ 


0.000000 


0.000000 


5 


0.000005 


1.1 X 10"'' 


0.000000 


0.000000 


6 


0.924624 


0.925554 


0.925550 


0.924957 


7 


0.000288 


0.000287 


0.000286 


0.000301 


8 


0.857835 


0.859467 


0.859457 


0.858464 


9 


0.000012 


1.4 X 10~^ 


0.000000 


0.000000 




U.UDDiZi 


U.UDDooD 


U.UDDo ( o 


U.UDD (94 


11 


0.074001 


0741 71 


0741 64 


n 074742 


12 


0.074001 




OTzLl fi4 

U.U / ^\.\J^ 


U.U / ^ / 


13 


0.000005 


1.1 X 10-^ 


0.000000 


0.000000 


14 


0.924624 


0.925554 


0.925550 


0.924957 


15 


0.008069 


0.008077 


0.008072 


0.008249 


16 


1.000000 


1.000000 


1.000000 


1.000000 








O22 




fit value 


1.704 


-1.408 


2.816 




evaluated value 


1.701 


-1.403 


2.806 



TABLE I. Comparison between ratios of determinants Ri and ratios of theta functions 9i for 
the genus two lattice with shape {rrii) = (11111). Ratios of determinants are given both for the 
largest Ay evaluation and for the Ay — 00 extrapolation. The theta function ratios are for a 
numerically fitted period matrix and for a period matrix evaluated according to the procedure of 
Section IV. 



lattice: 


cube 


L-shaped 


torus (++) 


genus 2 ( 1 — ) 


shape: 


(111) 


(121) 


(11111) 


(11) 


(11111) 


fits 


-0.19468 


0.19043 


-0.20531 


0.00080 


0.12858 




-0.19463 


0.19142 


-0.20573 


0.00062 


0.12786 




-0.19459 


0.19201 


-0.20571 


0.00050 


0.12735 




-0.19456 


0.19240 


-0.20567 


0.00041 


0.12698 


expected: 


-0.19444 


-0.20556 





0.125 



TABLE II. The logarithmic correction C for lattices with various topologies. A sequence of fits 
of four consecutive lattices of increasing size is seen to converge smoothly to the predicted value. 
For a more complete set of data see^. 
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# squares 


e 


6* 




3 


7r/2 


37r/2 


-0.024306 


5 


-7r/2 


57r/2 


0.01875 


8 


-27r 


47r 


0.0625 



TABLE III. The deficit angle, spanning angle and conical singularity contribution for vertices 
with various coordination numbers. 
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